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ABSTRACT 

We calculate the minimum mass of heavy elements required in the envelopes 
of Jupiter and Saturn to match the observed oversolar abundances of volatiles. 
Because the clathration efficiency remains unknown in the solar nebula, we have 
considered a set of sequences of ice formation in which the fraction of water 
available for clathration is varied between and 100%. In all the cases consid- 
ered, we assume that the water abundance remains homogeneous whatever the 
heliocentric distance in the nebula and directly derives from a gas phase of solar 
composition. Planetesimals then form in the feeding zones of Jupiter and Sat- 
urn from the agglomeration of clathrates and pure condensates in proportions 
fixed by the clathration efficiency A fraction of Kr and Xe may have been se- 
questrated by the Hjj" ion in the form of stable XeHjj" and KrHjj" complexes in 
the solar nebula gas phase, thus implying the formation of at least partly Xe- 
and Kr-impoverished planetesimals in the feeding zones of Jupiter and Saturn. 
These planetesimals were subsequently accreted and vaporized into the hydrogen 
envelopes of Jupiter and Saturn, thus engendering volatiles enrichments in their 
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atmospheres, with respect to hydrogen. Taking into account both refractory and 
volatile components, and assuming plausible molecular mixing ratios in the gas 
phase of the outer solar nebula, we show that it is possible to match the observed 
enrichments in Jupiter and Saturn, whatever the clathration efficiency Our cal- 
culations predict that the O/H enrichment decreases from ~ 6.7 to 5.6 times 
(O/H) in the envelope of Jupiter and from 18.1 to 15.4 times (O/H) in the 
envelope of Saturn with the growing clathration efficiency in the solar nebula. 
As a result, the minimum mass of ices needed to be injected in the envelope of 
Jupiter decreases from ~ 14.9 to 12.6 M®, including a mass of water diminishing 
from 9.6 to 7.7 M®. In the same conditions, the minimum mass of ices needed 
in the envelope of Saturn decreases from ~ 12.2 to 10.6 M®, including a mass of 
water diminishing from 7.7 to 6.4 M®. The accretion of planetesimals with ices 
to rocks ratios ~ 1 in the envelope of Jupiter, namely a value derived from the 
bulk densities of Ganymede and Callisto, remains compatible with the mass of 
heavy elements predicted by interior models. On the other hand, the accretion of 
planetesimals with similar ice-to-rock in the envelope of Saturn implies a mass 
of heavy elements greater than the one predicted by interior models, unless a 
substantial fraction of the accreted rock sedimented onto the core of the planet 
during its evolution. 

Subject headings: stars: planetary systems - stars: planetary systems: formation 
- solar system: formation 

1. Introduction 

Measurements by the mass spectrometer aboard the Galileo probe have shown that the 
abundances of C, N, S, Ar, Kr and Xe are all enriched by similar amounts with respect 
to their solar abundances in the atmosphere of Jupiter (Mahaffy et al. 2000; Wong et al. 
2004). Moreover, recent Cassini CIRS observations have also confirmed what was previously 
inferred from ground-based measurements, that C is substantially enriched in the atmosphere 
of Saturn (Flasar et al. 2005; Fletcher et al. 2008). 

In order to interpret these enrichments, it has been proposed by Gautier et al. (2001a,b), 
Alibert et al. (2005a,b), Mousis et al. (2006) and Hersant et al. (2004,2008) that the main 
volatile compounds initially existing in the solar nebula gas phase were essentially trapped 
by crystalline water ice in the form of clathrates or hydrates in the feeding zones of Jupiter 
and Saturn. These ices then agglomerated and formed planetesimals that were ultimately 
accreted by the forming Jupiter and Saturn. This is the fraction of these icy planetesimals 
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that vaporized when entering the envelopes of the two growing planets which engendered 
the observed volatiles enrichments. 

All the authors cited above postulate a full efficiency for the clathration process in 
the solar nebula, implying that guest molecules had the time to diffuse through porous 
water-ice planetesimals in the solar nebula. This remains plausible only if collisions between 
planetesimals have exposed essentially all the ice to the gas over time scales shorter or 
equal to planetesimals lifetimes in the nebula (Lunine & Stevenson 1985). However, the 
efficiency of collisions between planetesimals to expose all the "fresh" ice over such a time 
scale still remains questionable and we have no evidence that clathration was important in 
the primordial nebula (Owen et al. 1999,2000). Moreover, in all the afore-mentioned works, 
the abundance of available water ice is considered as a free parameter to select the amount of 
volatile species that are fully enclathrated and that contribute to the observed enrichments. 
As a result, this leads to abundances of available water in the feeding zones of Jupiter and 
Saturn which depart significantly from that predicted from a solar composition gas ,and the 
validity of such choices still needs to be investigated. 

Here, we show that it is possible to explain the volatiles enrichments in Jupiter and 
Saturn by postulating an incomplete clathration process in the solar nebula, and also by using 
an abundance of available water derived from a homogeneous gas phase of solar composition. 
In the extreme case of no clathration in the outer nebula, it is still possible to match the 
enrichments observed in the two planets. Volatile species that are not enclathrated essentially 
form pure condensates at lower temperatures in the feeding zones of Jupiter and Saturn and 
icy planetesimals ultimately accreted by the two forming planets then agglomerate from 
clathrates and pure condensates in proportions fixed by the clathration efficiency. We finally 
show that the minimum amounts of ices needed in the envelopes of Jupiter and Saturn 
to match their observed enrichments are compatible with the amounts of heavy elements 
predicted by interior models. 

Sec. |2]is devoted to the description of the formation mechanisms of Jupiter and Saturn 
in the framework of the core-accretion model. We also briefly review the distribution of 
heavy elements in the two planets which is predicted by interior models. In Sec. [31 we 
describe the formation sequence of icy planetesimals in the feedings zones of Jupiter and 
Saturn as a function of the clathration efficiency. In Sec. HJ we determine the minimum 
masses of heavy elements required to match the oversolar abundances of volatiles observed 
in Jupiter and Saturn. Sec. [5] is devoted to discussion. 
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2. Delivery of volatiles to the envelopes of Jupiter and Saturn 

There are two main planet formation models considered in the recent literature. The 
first one, the disk instability model, postulates that planets form by gravitational instability 
in a protoplanetary disk. In the second model, the core-accretion model, planets form in 
two succesive phases. A solid core is first formed by collisional accretion of planetesimals, in 
a process similar to the formation of terrestrial planets. In a second phase, when the mass 
of this core becomes large enough, rapid gas accretion is triggered, leading to the formation 
of a gas giant planet. Recent improvements in core-accretion models (see e.g. Alibert et al. 
2005c; Hubickyj et al. 2005) allow the formation of gas giant on time scales compatibles 
with disk observed lifetimes. Moreover, by including migration and disk evolution, some of 
these models allow forming planets whose internal structures are comparable with the ones 
of Jupiter and Saturn (Alibert et al. 2005b). Finally, in a population synthesis approach, 
the same model deliver planet statistical properties consitent with observed ones (at least 
for single planets in quasi-circular orbits around G stars; see e.g. Mordasini et al. 2008). 

Our present work is based on the afore-mentioned extended core accretion model, taking 
into account migration and disk evolution. In this model, the two planets start their for- 
mation at larger distance from the Sun, migrate and simultaneously accrete gas and solids. 
When the disk has disappeared, the accretion stops, and the two planet do not migrate 
anymore. In this model, the structure of the protoplanetary disk is derived assuming that 
viscous heating is the predominant heating source. Under this hypothesis, disk models show 
that the outer parts of the disk are protected from solar irradiation by shadowing effect 
of the inner disk parts, and temperature in the planet forming region (between 5 and 15 
AU) decreases down to very low values (~ 20 K). However, note that irradiation onto the 
central disk parts could modify the disk structure so much that shadowing effect would be 
prevented in the outer parts. In this case, the temperature in the planet-forming region 
would be higher. 

In this paper, we compare the amount of heavy elements needed to explain the enrich- 
ments in volatiles species in Jupiter and Saturn with theoretical determination of the planets 
metal content. In this context, we use the internal structure models of Saumon & Guillot 
(2004) (hereafter SG04) who have derived estimates for the core mass and mean metallicity 
of Jupiter and Saturn, the main uncertainty in these determinations being the equation of 
state of hydrogen and helium under high pressure. Taking into account these different un- 
certainties, the total mass of heavy elements present in Jupiter (M corc + Mz iCnvc ) can be as 
high as ~ 42 M e whereas the mass of the core can range between and 13 M®. In the case 
of Saturn, the mass of heavy elements can increase up to ~ 35 M @ with the envelope mass 
varying between ~ and 10 M e and the core mass ranging between ~ 8 and 25 M e . Note 
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that these estimates will be improved, in the case of Jupiter, by the future JUNO mission. 

3. Formation sequence of icy planetesimals in the giant planets feeding zones 

We describe here the formation sequence of the different ices produced in the feeding 
zones of proto- Jupiter and proto-Saturn. Once formed, these ices will add to the composition 
of the planetesimals that will be accreted by the two giant planets during their growth. 

In order to define the initial gas phase composition of the solar nebula, we assume 
that the abundances of all elements, including oxygen, are protosolar (Lodders 2003) and 
consider both refractory and volatile components. Refractory components include rocks 
and organics. According to Lodders (2003), rocks contain ~ 23% of the total oxygen in 
the nebula. The fractional abundance of organic carbon is assumed to be 55% of total 
carbon (Pollack et al. 1994), and the ratios of C/O/N included in organics is assumed to be 
1/0.5/0.12 (Jessberger et al. 1988). We then assume that the remaining O, C, and N exist 
only under the form of H 2 0, CO, CO2, CH 3 OH, CH 4 , N 2 , and NH 3 . Hence, once the gas 
phase abundances of elements are defined, the abundances of CO, CO2, CH 3 OH, CH 4 , N 2 
and NH 3 are determined from the adopted CO/CO2/CH3OH/CH4 and N 2 /NH 3 gas phase 
molecular ratios, and from the C/O/N relative abundances set in organics. Finally, once 
the abundances of these molecules are fixed, the remaining O gives the abundance of H 2 0. 
We then set CO/C0 2 /CH 3 OH/CH 4 = 70/10/2/1 in the gas phase of the disk, values that 
are consistent with the ISM measurements considering the contributions of both gas and 
solid phases in the lines of sight (Frerking et al. 1982; Ohishi et al. 1992; Ehrenfreund & 
Schutte 2000; Gibb et al. 2000). In addition, S is assumed to exist in the form of H 2 S, 
with n 2 S/n 2 = 0.5 x (S/H 2 ) Q , and other refractory sulfide components (Pasek et al. 2005). 
We also consider N 2 /NH_3 = 1/1 in the nebula gas-phase. This value is compatible with 
thermochemical calculations in the solar nebula that take into account catalytic effects of 
Fe grains on the kinetics of N 2 to NH 3 conversion (Fegley 2000). In the following, we adopt 
these mixing ratios as our nominal model of the solar nebula gas phase composition (see 
Tabled]). 

The process by which volatiles are trapped in icy planetesimals, illustrated in Figs. [1] 
and [2j is calculated using the stability curves of hydrates, clathrates and pure condensates, 
and the thermodynamic path detailing the evolution of temperature and pressure at 5.2 and 
9.5 AU in the solar nebula, corresponding to the actual positions of Jupiter and Saturn, 
respectively. We refer the reader to the works of Papaloizou & Terquem (1999) and Alibert 
et al. (2005c) for a full description of the turbulent model of accretion disk used here. 
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The stability curves of hydrates and clathrates derive from Limine Sz Stevenson (1985) 's 
compilation of published experimental work, in which data are available at relatively low 
temperatures and pressures. On the other hand, the stability curves of pure condensates 
used in our calculations derive from the compilation of laboratory data given in the CRC 
Handbook of Chemistry and Physics (Lide 2002). The cooling curve intercepts the stability 
curves of the different ices at particular temperatures and pressures. For each ice considered, 
the domain of stability is the region located below its corresponding stability curve. The 
clathration process stops when no more crystalline water ice is available to trap the volatile 
species. Note that, in the pressure conditions of the solar nebula, CO2 is the only species 
that crystallizes at a higher temperature than its associated clathrate. We then assume that 
solid CO2 is the only existing condensed form of CO2 in this environment. In addition, 
we have considered only the formation of pure ice of CH 3 OH in our calculations since, to 
our best knowledge, no experimental data concerning the stability curve of its associated 
clathrate have been reported in the literature. 

Because the clathration efficiency remains unknown in the solar nebula, we have con- 
sidered a set of sequences of ices formation in which the fraction of water available for 
clathration is varied between and 100%. Figure [T] illustrates the case where the efficiency 
of clathration is total, implying that guest molecules had the time to diffuse through porous 
water-ice planetesimals before their accretion by proto- Jupiter and proto-Saturn. In this 
case, NH 3 , H 2 S, PH 3 , Xe, CH 4 and ~ 61% of CO form NH 3 -H 2 hydrate and H 2 S- 
5.75H 2 0, PH 3 -5.67H 2 0, Xe-5.75H 2 0, CH 4 -5.75H 2 and CO-5.75H 2 clathrates with the 
available water. The remaining CO, as well as N 2 , Kr, and Ar, whose clathration normally 
occurs at lower temperatures, remain in the gas phase until the nebula cools enough to al- 
low the formation of pure condensates. Figure [2] illustrates the case where the efficiency of 
clathration is only ~ 10%. Here, either only a part of the clathrates cages have been filled 
by guest molecules, either the diffusion of clathrated layers through the planetesimals was 
to slow to enclathrate most of the ice, or the poor trapping efficiency was the combination 
of these two processes. In this case, only NH 3 and ~ 21% of H 2 S form NH 3 -H 2 hydrate 
and H 2 S-5.75H 2 clathrate. Due to the deficiency in accessible water in icy planetesimals, 
the remaining H 2 S and all PH 3 , Xe, Kr, CH 4 , CO, Ar and N 2 form pure condensates in the 
solar nebula. 

Table [2] summarizes the trapping/formation conditions of the different ices 
calculated at 5.2 and 9.5 AU in the solar nebula in the cases of 100% and 10% 
clathration efficiencies, and for our nominal gas phase. Using these thermody- 
namic conditions, one can estimate the mass ratios of the different ices with 
respect to H 2 in the planetesimals formed in the solar nebula. Indeed, the 
volatile, i, to water mass ratio in these planetesimals is determined by the rela- 
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tion given by Mousis & Gautier (2004): 



rrii = 



X, E(r;Ti,Pi) 



(1) 



Xh 2 o ^(^; T H2 o, Ph 2 o) 



where and Xh 2 o are the mass mixing ratios of the volatile i and H 2 with 
respect to H 2 in the nebula, respectively. S(i?;T i5 Pj) and H(R;T H2 o, Ph 2 o) are 
the surface density of the nebula at a distance r from the Sun at the epoch of 
hydratation or clathration of the species i, and at the epoch of condensation 
of water, respectively. From m i5 it is possible to determine the mass fraction 
Mi of species i with respect to all the other volatile species taking part to the 
formation of an icy planetesimal via the following relation: 



Note that, whatever their formation distance in the nebula, the composition of plan- 
etesimals remains almost constant, provided that the gas phase composition does not vary 
(Marboeuf et al. 2008) and that the clathration efficiency remains constant. In 
particular, if they formed in the same gas phase conditions and with the same 
clathration efficiency in the nebula, the ices accreted by proto-Jupiter and proto- 
Saturn share a similar composition along their migration pathways. On the other 
hand, the composition of ices can be somewhat altered by the variation of the 
clathration efficiency (see Table [3]). This does not impair the quality of the fits 
of the observed enrichments but the budget of heavy elements needed in the 
envelopes of the two planets can be modified (see Sec. 3]). However, in any case, 
it appears difficult to provide some constraint on the migration status of Jupiter 
and Saturn from their observed enrichments because the changes of planetesi- 
mals composition are never significant in the giant planets formation zone, even 
when the clathration efficiency is varied. 



Mi 




(2) 



j=l,n 



with 
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4. Minimum masses of heavy elements in the envelopes of Jupiter and Saturn 

From the adjustment of the masses of ices vaporized in the envelopes, we have been 
able to reproduce the observed volatiles enrichments. Our strategy was i) to match 
the maximum number of observed volatiles enrichments and ii) in the range of 
possible solutions, to select the fit that minimizes the mass of ices needed in the 
envelopes of Jupiter and Saturn. 

In our calculations, we have considered the measurements of C, N, S, Ar, 
Kr and Xe abundances in Jupiter's atmosphere determined using the mass spec- 
trometer on board the Galileo probe (Mahaffy et al. 2000; Wong et al. 2004). 
Estimates of the abundance of P from 10 /im Cassini/CIRS PH 3 observations 
have been recently updated by Fletcher et al. (2009) to include improved esti- 
mates of the mid-IR aerosol opacity, spectroscopic linedata and NH3 abundance, 
resulting in derived abundances slightly larger than the previous analysis of Irwin 
et al. (2004). 

In the case of Saturn, we have only used the recent determination of C 
abundance by Fletcher et al. (2008) from Cassini CIRS spectra. This new mea- 
surement updates the previous study by Flasar et al. (2005) using thousands 
of high-resolution spectra in both the mid- and far-infrared acquired during 
Cassini's prime mission, taken at a range of spatially-resolved locations on the 
planet. As well as improving the precision of the measurement, Fletcher et al. 
(2008) demonstrate the consistency between results obtained from rotational and 
vibrational line manifolds in the far- and mid-IR, and show the lack of hemispher- 
ical asymmetry in the CH 4 abundance, confirming the hypothesis that this gas 
is well-mixed throughout the observable atmosphere and is therefore represen- 
tative of the bulk composition. Finally, we use a new estimate of Saturn's PH3 
abundance, taking into account the latitudinal variability and mid-IR aerosol 
opacities derived by Fletcher et al. (2009). 

We have chosen to omit the mixing ratios of condensible species (NH 3 , H 2 S), 
as infrared remote sensing is unable to constrain the abundances beneath the 
condensation clouds, and the competing spectral effects of these species are 
difficult to disentangle from measurements of Saturn's radio frequency opacity 
(Briggs & Sackett 1989). Furthermore, the absence of an atmospheric entry- 
probe at Saturn means we have no constraints on the abundances of the noble 
gases (Ar, Kr, Xe) in that planet's atmosphere. 

Table H] gives our "minimum" fits for the two planets in the case of the nominal model 
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and with 100% and 10% clathration efficiencies. In the case of Saturn and whatever 
the clathration efficiency, the observed C and P enrichments are systematically 
matched within error bars or range of values. Similarly, in the case of Jupiter, 
the observed C, N, S and Ar enrichments are systematically matched within 
error bars. However, the calculated P abundance in the Jovian atmosphere is 
higher than the observed value. On the other hand, the relationship between 
the PH 3 abundance observed in the 1—4 bar region on both planets and the 
P/H content of the interior is dependent upon the assumed value of the eddy 
mixing coefficient at the kbar level. It is therefore possible that the observed 
PH 3 abundances provide only lower bounds on the P/H abundance (Fegley & 
Lodders 1994; Fletcher et al. 2009). Moreover, the calculated Kr and Xe enrichments 
are higher than the observed ones but the presence of H3 ion in the primitive nebula, which 
induces an efficient trapping of these species in the form of stable complexes XH^ (with 
X = Kr and Xe; Pauzat Sz Ellinger 2007) in the gas phase, may limit their ability to be 
enclathrated or to condense in the outer nebula (Mousis et al. 2008). Therefore, if the 
abundance were comparable to those of Kr and Xe in Jupiter's feeding zone (H3 /H2 ~ 
10 -9 -10~ 10 ), which is a reasonable assumption (see e.g. Mousis et al. 2008), the resulting 
enrichments of these two nobles gases in the Jovian envelope should be lower than the values 
calculated here and might match the observed values^. For example, if H^/H^ ~ 5 x 10~ 10 
in the feeding zone of Jupiter and if similar amounts of Kr and Xe have been trapped by , 
their corresponding enrichments are now 2.4 and 1.7 and match the observed values in the 
case of full clathration efficiency. Note that, since the H3" abundance is expected to increase 
with the growing heliocentric distance, the observed deficiency of Titan's atmosphere in 
Kr and Xe (Niemann et al. 2005) was suggested to be caused by the presence of a higher 
concentration of KrHjj" and XeH^ complexes in Saturn's feeding zone, inducing the formation 
of Kr- and Xe-poor planetesimals ultimately accreted by the satellite (Mousis et al. 2008). If 
this scenario is correct, the abundances of Kr and Xe should be solar in Saturn's atmosphere. 
In a less extreme case, if we also adopt /H2 ~5 x 10~ 10 in Saturn's feeding zone, Kr and 
Xe enrichments should be of 6.8 and 4.7, instead of respectively 7.4 and 11.6 in the case of 
full clathration efficiency. 

Independently of the efficiency of Xe and Kr trapping by H3 in the nebula^, our calcu- 



1 Ar also forms stable ArHj complexe in the gas phase (Pauzat & Ellinger 2007) but its abundance in 
the solar nebula is several orders-of-magnitude higher than the one envisaged for H J . Hence, the fraction of 
Ar sequestrated by can be neglected. 

2 Kr and Xe poorly influence the budget of ices accreted by Jupiter and Saturn, due to their low abundances 
compared to those of other volatiles. 
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lations predict that the O/H enrichment decreases from ~ 6.7 to 5.6 times (O/H) in the 
envelope of Jupiter and from ~ 18.1 to 15.4 times (O/H) in the envelope of Saturn, with the 
growing clathration efficiency in the solar nebula. Figure [3] shows that the minimum masses 
of ices required in the envelopes of Jupiter and Saturn to match the observed enrichments 
diminish when the clathration efficiency increases in the solar nebula. This effect illustrates 
that the variation of the clathration efficiency affects the composition of ices produced in 
the outer nebula. The minimum mass of ices thus needed to be injected in the envelope of 
Jupiter decreases from ~ 14.9 to 12.6 M®, including a mass of water that diminishes from 
9.6 to 7.7 M®, with the growing clathration efficiency. In the same conditions, the minimum 
mass of ices needed in the envelope of Saturn decreases from ~ 12.2 to 10.6 M®, including 
a mass of water that diminishes from 7.7 to 6.4 M®. 

5. Summary and discussion 

• In this report, we have assumed that the water abundance is homogeneous whatever 
the considered heliocentric distance and derives directly from a gas phase of solar 
composition. The amount of water available for clathration is then no more a free 
parameter adjusted in order to ease the fits of the observed volatiles enrichments in 
Jupiter and Saturn; 

• Planetesimals then form in the feeding zones of Jupiter and Saturn from clathrates 
and pure condensates in proportions fixed by the clathration efficiency; 

• A fraction of Kr and Xe may have been sequestrated by in the form of XeHjj" and 
KrH^" complexes in the solar nebula gas phase, thus implying the formation of at least 
partly Xe- and Kr-depleted planetesimals in the feeding zones of Jupiter and Saturn; 

• From plausible molecular mixing ratios in the gas phase of the outer solar nebula and 
from the calculation of the composition range of ices accreted by the growing Jupiter 
and Saturn, we show that it is possible to match the observed C, N, S, Ar, Kr and 
Xe enrichments, whatever the clathration efficiency. 

The minimum mass of ices [12.6-14.9 M®] required in the envelope of Jupiter to match 
the observed enrichments is lower than the maximum amount of heavy elements (42 M®) 
predicted in the same zone by the internal structure models of SG04. Even if planetesimals 
accreted by proto- Jupiter are half composed of ices and half made of rocks (ices to rocks ratio 
of 1), a value compatible with the internal structures of Ganymede and Callisto (Sohl et al. 
2002), the global mass of solids injected in the envelope remains lower than the maximum 
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one predicted by SG04. These results are then compatible with formation scenarios of the 
two Galilean satellites from the accretion of planetesimals formed in the nebula without 
having been vaporized inside the subdisk (Mousis & Gautier 2004; Mousis & Alibert 2006). 

On the other hand, the minimum mass of ices [10.6-12.2 M®] needed in the envelope of 
Saturn to match the observed enrichments slightly exceeds the maximum amount of heavy 
elements (~ 10 M®) predicted by SG04. However, note that the minimum mass of ices 
predicted by our model is lower than the [13.7-18.3M®] range of values derived by Mousis et 
al. (2006) from the assumption that crystalline water was abundant enough in the feeding 
zones of Jupiter and Saturn to enclathrate essentially all the volatile species present in the 
gas phase. 

It is still possible to slightly decrease the minimum mass of ices needed in the envelopes 
of Jupiter and Saturn by adopting C0 2 /CO gas phase ratios different from the one selected in 
our nominal model. The variation of this ratio rules the distribution of C and O among CO2, 
CO, CH3OH, CH 4 and H 2 molecules and then strongly affects the budget of crystalline 
water available for clathration. Figure H] illustrates this effect and shows that the minimum 
mass of ices needed in the envelopes of Jupiter and Saturn to match the observed enrichments 
diminishes when the CO2/CO ratio increases in the initial gas phase of disk. In the case of 
Saturn, CO2/CO ratios greater than ~ 0.2 in the nebula gas phase allow the minimum mass 
of ices needed in the envelope to be lower than the maximum value given by SG04. However, 
it is difficult to explain the accretion in Saturn of planetesimals with ice to rock ratios ~ 
1 -similar to that predicted in Titan by internal structure models (Tobie et al. 2006)- by 
invoking higher CO2/CO ratios in the nebulae gas phase, in agreement with the amount of 
ices predicted in the envelope by SG04. 

Indeed, the extremes adopted in Fig. H]for the C0 2 /CO ratio in the nebula -0.1-1.0- 
correspond to contributions from, respectively, the ISM gas and solid phases (Ehrenfreund 
& Schutte 2000; Gibb et al. 2004). Hence, any CO2/CO mixing ratio adopted in the solar 
nebula gas phase should hold within this range of values. On the other hand, if a substantial 
fraction of rock contained in planetesimals that accreted in the envelope has sedimented 
onto the core of Saturn during its evolutional, the inconsistency between the accreted mass 
of planetesimals and the one predicted by interior models could be removed. 

Here, we have assumed that the observed enrichments in volatiles were en- 
gendered by the vaporization of icy planetesimals when they entered the en- 
velopes of the growing planets. On the other hand, the erosion of a significant 



3 Due to their higher density and their low volatility, rock should reach the deepest layers of Saturn's 
envelope during the accretion of planetesimals and remain in the solid phase at extremely high pressure. 
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part of the giant planets cores could also constitute a possible source of volatiles 
enrichments in their atmospheres. However, in this case, the solids accreted by 
the cores should share the same composition as those accreted later in the en- 
velope. Hence, the approach described in this work remains valid whatever the 
delivery process of volatiles (core erosion or accretion of solids in the envelope). 

Finally, we note that a new compilation of protosolar abundances has been recently 
published by Grevesse et al. (2007). Comparisons between this compilation and the one 
used in the present work and taken from Lodders (2003) show that, except for Ar abundance 
which is ~ 2.5 times smaller in the more recent tabulation, no such substantial variation 
is observed. If Ar is not considered, the use of the compilation of Grevesse et 
al. (2007) do not alter the conclusions of this paper. However, if we consider the 
Ar abundance given by Grevesse et al. (2007), only the observed C, N and S enrichments 
in Jupiter can be matched by our calculations, whatever the clathration efficiency The 
use of two different methods seems to be at the origin of the difference in the quoted solar 
Ar abundance in these compilations. Indeed, Lodders (2003) based the Ar abundance on 
nucleosynthesis arguments whereas Grevesse et al. (2007) based the Ar abundance from 
abundance ratios of Ar to other elements measured in the solar wind (SW) and solar energetic 
particles (SEPs) coming from the Sun's corona. Following Lodders (2008), the limitation of 
the method employed by Grevesse et al. (2007) is that elements may become fractionated in 
the SW and SEPs from photospheric abundances according to their different first ionization 
potentials (FIP). Relative to photospheric values, elements with low FIP < 10 eV would be 
enriched in the SW and SEPs, whereas elements with high FIP (such as Ar) appear to be 
depleted (Lodders 2003,2008). From these considerations, and for reasons of consistency, we 
have adopted the compilation of Lodders (2003) in all our calculations. 
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Table 1: Elemental and molecular abundances in the solar nebula 



Species X 


X/H 2 




C 


5.82 x 10" 


-4 


N 


1.60 x 10" 


-4 





1.16 x 10" 


-3 


s 


3.66 x 10" 


-5 


Ar 


8.43 x 10" 


-6 


Kr 


4.54 x 10" 


-9 


Xe 


4.44 x 10- 


10 


P 


6.88 x 10" 


-7 


H 2 


4.44 x 10" 


-4 


CO 


2.21 x 10" 


-4 


co 2 


3.16 x 10" 


-5 


NH 3 


4.05 x 10" 


-5 


H 2 S 


1.83 x 10" 


-5 


N 2 


4.05 x 10" 


-5 


CH 3 OH 


6.31 x 10" 


-6 


CH 4 


3.16 x 10" 


-6 


PH 3 


6.88 x 10" 


-7 


Kr 


4.54 x 10" 


-9 


Xe 


4.44 x 10- 


10 



Note. — Elemental abundances derive from Lodders (2003). Molecular abundances are determined from 
our nominal gas phase composition. 
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Table 2: Thermodynamic conditions of ices formation in the outer nebula 



100% clathration efficiency 


Heliocentric distance (AU) 




c 
o.z 








9.5 






Ice 


T(K) 


P (bars) 


E (g.cm 2 ) 


T(K) 


P (bars) 




E (g.cm 2 ) 


H9O 


155.65 


3 35 x 10" 7 


670.73 


155.51 


2.03 


x 10- 


-7 


1146.91 


CH3OH 


103.18 


2.18 x 10~ 7 


536.40 


104.35 


1.32 


x 10" 


-7 


915.31 


NH 3 -H 2 


88.31 


1.87 x 10~ 7 


495.52 


88.65 


1.13 


x 10" 


-7 


845.97 


H 2 S-5.75H 2 


83.89 


1.77 x lO" 7 


482.24 


85.23 


1.08 


x 10- 


-7 


829.15 


co 2 


74.55 


1.57 x 10~ 7 


452.59 


75.76 


9.62 


x 10- 


-8 


780.15 


PH 3 -5.67H 2 


70.84 


1.49 x 10~ 7 


440.10 


70.05 


8.89 


x 10- 


■8 


748.65 


Xe-5.75H 2 


59.83 


1.25 x 10- 7 


400.05 


59.95 


7.58 


x 10- 


-8 


688.23 


CH 4 -5.75H 2 


54.37 


1.13 x 10~ 7 


378.12 


55.49 


6.99 


x 10- 


-8 


659.24 


CO-5.75H 2 


47.46 


9.72 x 10~ 8 


347.84 


47.59 


5.95 


x 10- 


-8 


603.46 


Kr 


29.15 


5.43 x 10~ 8 


244.10 


29.46 


3.45 


x 10- 


-8 


439.45 


CO 


25.33 


4.47 x 10~ 8 


214.44 


25.67 


2.90 


x 10- 


-8 


393.47 


Ar 


22.17 


3.64 x 10~ 8 


185.64 


22.53 


2.42 


x 10- 


-8 


349.41 


N 2 


21.60 


3.49 x 10- 8 


179.95 


21.83 


2.32 


x 10- 


-8 


338.70 


10% clathration efficiency 


Heliocentric distance (AU) 




5.2 








9.5 






Ice 


T(K) 


P (bars) 


E (g.cm" 2 ) 


T(K) 


P (bars) 




E (g.cm" 2 ) 


H 2 


155.65 


3.35 x 10- 7 


670.73 


155.51 


2.03 


x 10- 


-7 


1146.91 


CH3OH 


103.18 


2.18 x 10~ 7 


536.40 


104.35 


1.32 


x 10- 


-7 


915.31 


NH 3 -H 2 


88.31 


1.87 x 10- 7 


495.52 


88.65 


1.13 


x 10- 


-7 


845.97 


H 2 S-5.75H 2 


83.89 


1.77 x 10~ 7 


482.24 


85.23 


1.08 


x 10- 


-7 


829.15 


co 2 


74.55 


1.57 x 10~ 7 


452.59 


75.76 


9.62 


x 10- 


-8 


780.15 


H 2 S 


68.80 


1.45 x 10~ 7 


433.01 


70.05 


8.89 


x 10- 


-8 


748.65 


PH 3 


45.52 


9.28 x 10~ 8 


338.73 


45.81 


5.71 


x 10- 


-8 


589.95 


Xe 


38.20 


7.59 x 10~ 8 


301.06 


37.96 


4.64 


x 10- 


-8 


524.77 


Kr 


29.15 


5.43 x 10~ 8 


244.10 


29.46 


3.45 


x 10- 


-8 


439.45 


CH 4 


27.87 


5.11 x 10~ 8 


234.69 


28.44 


3.30 


x 10- 


-8 


427.78 


CO 


25.33 


4.47 x 10~ 8 


214.44 


25.67 


2.90 


x 10- 


-8 


393.47 


Ar 


22.17 


3.64 x 10~ 8 


185.64 


22.53 


2.42 


x 10- 


-8 


349.41 


N 2 


21.60 


3.49 x 10" 8 


179.95 


21.83 


2.32 


x 10" 


-8 


338.70 



Note. — T, P and S are the temperature, pressure and surface density of the H 2 -dominated gas given 
at 5.2 or 9.5 AU in the solar nebula at the epoch of condensation of volatile i or its trapping by crystalline 



Table 3: Average composition of ices formed in the outer solar nebula 



Ice 


100% clathration efficiency 


10% clathration efficiency 


H 2 


6.08 x 10- 1 


6.39 x 10- 1 


CO 


2.00 x 10- 1 


1.64 x 10- 1 


co 2 


7.18 x 10~ 2 


7.54 x 10~ 2 


NH 3 


3.88 x 10" 2 


4.07 x 10~ 2 


H 2 S 


3.42 x 10~ 2 


3.30 x 10~ 2 


N 2 


2.44 x 10- 2 


2.56 x 10~ 2 


CH3OH 


1.23 x 10~ 2 


1.29 x 10~ 2 


Ar 


6.73 x 10~ 3 


7.07 x 10~ 3 


CH 4 


2.19 x 10~ 3 


1.46 x 10- 3 


PH 3 


1.17 x 10~ 3 


9.55 x 10~ 4 


Kr 


1.08 x 10~ 5 


1.14 x 10~ 5 


Xe 


2.66 x 10~ 6 


2.11 x 10~ 6 



Note. — Composition of ices is calculated in the cases of 100% and 10% clathration efficiencies. Ratio of 
the mass of ice i to the global mass of ices (wt%) is determined from our nominal gas phase composition. 
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Table 4: Observed enrichments in volatiles in Jupiter and Saturn, and calculated enrichments 
in the case of our nominal model 





Jupiter 






Saturn 






Species 


Observed 


(1) 


(2) 


Observed 


(1) 


(2) 







5.6 


6.7 




15.4 


18.1 


c 


4.1 ± r 


3.1 


3.1 


9.2±0.4 d 


8.8 


8.8 


N 


4.15 ± 1.6° 


3.0 


3.7 




8.6 


10.4 


S 


2.4 ±0.6" 


2.5 


2.9 




7.0 


7.8 


P 


3.2±0.15 6 


4.6 


4.4 


8.9 - 13.5 6 


12.6 


12.1 


Ar 


2.15±0.4 C 


1.9 


2.4 




5.9 


7.1 


Kr 


2 ± 0.4 C 


2.6 


3.2 




7.4 


9.0 


Xe 


2 ± 0.4 C 


4.2 


3.9 




11.6 


10.7 



Note. - Cases (1) and (2) correspond to 100% and 10% clathration efficiencies, respectively. The 
sequestration of Kr and Xc in the forms of KrH^~ and XeH^" in the solar nebula gas phase is not taken into 
account in the presented calculations (see text). The observed values are derived from a Wong et al. 2004, b 
Fletcher et al. (2009), c Mahaffy et al. (2000) and d Fletcher et al. (2008), using the protosolar abundances 
of Loddcrs (2003). 
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Fig. 1.— Stability curves of hydrate (NH 3 -H 2 0), clathrates (X-5.75H 2 or X-5.67H 2 0) (solid 
lines), and pure condensates (dotted lines), and cooling curves of the Solar nebula at the 
heliocentric distances of 5.2 and 9.5 AU, respectively, assuming a full efficiency of clathration. 
Abundances of various elements are solar, with CO/CO2/CH3OH/CH4 = 70/10/2/1, H 2 S/H 2 
= 0.5 x (S/H 2 ) , and N 2 /NH 3 = 1/1 in the gas phase of the disk. Species remain in the gas 
phase above the stability curves. Below, they are trapped as clathrates or simply condense. 
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Total pressure (bars) 



Fig. 2. — Same as in Fig. 1, but with a 10% clathration efficiency. 
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Fig. 3. — Minimum masses of ices and water needed in the envelopes of Jupiter and Saturn 
to fit the measured abundances of volatiles in the case of our nominal model, as a function 
of the clathration efficiency in the two giant planets feeding zones. 
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0.2 0.4 0.6 0.8 1 

C0 2 /CO 

Fig. 4. — Minimum masses of ices and water needed in the envelopes of Jupiter and Saturn 
to fit the measured abundances of volatiles, as a function of the CO2/CO ratio postulated in 
the solar nebula gas phase. Abundances of various elements are solar, with CO/CH3OH/CH4 
= 70/2/1, H 2 S/H 2 = 0.5 x (S/H 2 ) , and N 2 /NH 3 = 1/1 in the gas phase of the disk. The 
vertical arrow indicates the C0 2 /CO ratio corresponding to our nominal model. Calculations 
were made with the assumption of a full clathration efficiency. 



